Electronic and dynamic properties of a one-dimensional Thue-Morse chain 
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We study the electronic and dynamic properties of a one-dimensional Thue-Morse chain within 
the framework of the transfer model. By means of direct diagonalization of the Hamiltonian matrix 
we first show that the trace map of the transfer matrices is exactly the same as that in the diagonal 
model. Then, by calculating several quantities such as the wave function, the Lyapunov exponent, 
and the Landauer resistivity, we show that all the electronic states are extended despite the singular 
continuity of the energy spectrum. Our results indicate that the electronic properties of the Thue- 
Morse chain is independent of the kind of the model, which is contrary to the result of Chakrabarti 
et al. [Phys. Rev. Lett. 74, 1403 (1995)]. To deepen our understanding, we study the dynamics 
of an electronic wave packet and show that the wave packet spreads superdiffusively for long times 
with dynamic indices intermediate between periodic and quasiperiodic chains. Several features of 
dynamics distinctive from the case of the Fibonacci chain are discussed. We also study the effects 
of electron-phonon interaction on the dynamics of the wave packet by taking into account a kind 
of nonlinear interaction. Degree of dynamic localization is shown to be crucially dependent on the 
strength of the hopping energy and the nonlinear interaction parameter. 



PACS numbers: 71.25.-s, 61.44.+p, 64.60.Ak, 71. 20. Ad 



I. INTRODUCTION 

Since the discovery of the quasicrystalline phase 1 much 
attention has been devoted to the quasiperiodic sys- 
tems. The Fibonacci chain, a one-dimensional version 
of quasicrystal characterized by either the inflation rule 
(L,S) — ► (LS,L) or the substitution rule Si+i = SiSi-i 
with Si = {L}, is the most intensively studied example, 
and new concepts on the electronic properties such as 
the singular continuous energy spectrum and the critical 
states are now well established. 2 Stimulated by unusual 
electronic properties of the Fibonacci chain as well as the 
development of layer-growth technique, 3 there has been a 
variety of interest in one-dimensional deterministic aperi- 
odic chains. Among them the most well-known example 
is the Thue-Morse chain characterized by either the in- 
flation rule (L, S) — > (LS, SL) or the substitution rule 
Si+i = SiSi with Si = {£}. 

Unusual electronic and Fourier spectral properties of 
the Thue-Morse chain are often compared with those of 
the Fibonacci chain. According to the nature of the 
Fourier spectrum of the sequence, the Thue-Morse se- 
quence is more disordered than the Fibonacci sequence 
since the former has a singular continuous spectrum 4 
while the latter has a point-like spectrum with Bragg 
5-peaks. In the meanwhile, according to the nature of 
the electronic energy spectrum, both of the chains are 
located in the same category; both have the singular 
continuous energy spectrum, which is a general feature 
of one-dimensional deterministic aperiodic chains. 5 How- 
ever, according to the nature of the electronic states, the 
Thue-Morse chain is much closer to the periodic chain 
than the Fibonacci chain; all the electronic states of the 
latter are rigorously proven to be critical, 2 while most 
of the electronic states of the former exhibit Bloch-likc 



extended character. 6 

In studying the electronic properties 
dimensional deterministic aperiodic chains, it 
common to use the tight-binding equation 
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Here v n is the on-site energy of the site n, t n ±i,n is the 
nearest neighboring hopping energy between the sites n 
and nil, and T n is the transfer matrix. 

When a sequence to be studied is introduced into v n 
(tnm) with setting t nm — 1 (v n — 0), it is called the diag- 
onal (transfer) model. In the case of the Fibonacci chain, 
there is no preference in choosing the kind of the model 
and both of the models have been equally well studied. 
However, in the case of the Thue-Morse chain, most of 
the existing literatures have been devoted to the diag- 
onal model rather than the transfer model because the 
transfer matrix Mi = Il^ =1 T n appeared in the transfer 
model is so complex that the trace map method which is 
the most powerful tool in studying quasiperiodic or ape- 
riodic chains cannot be easily applied. 7 Furthermore, the 
complexity of the transfer matrix in the transfer model 
also leads a literature that contains incorrect argument 
on the electronic properties of the Thue-Morse chain: 
Chakrabarti et al. s argued that the transfer model does 
not support the extended states unlike to the case of the 
diagonal model. However, we argue that the electronic 
properties of the Thue-Morse chain are independent of 
the kind of the model. To certify our argument, we show 
that the trace map of the transfer matrices in the transfer 
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model is exactly the same as that in the diagonal model 
and that all the electronic states except the two edge 
states are extended. 

It has been generally believed that the electronic states 
corresponding to a singular continuous energy spectrum 
are critical. Of course, in some cases such as the gener- 
alized Fibonacci chain 9 and the period-doubling chain, 10 
existence of extended states has been reported. How- 
ever, even in these cases, the portion of the extended 
states to the critical states goes to zero in the infinite 
limit of the chain size, and the global electronic proper- 
ties of the chains are qualitatively similar to those of the 
Fibonacci chain. In this context, the Thue-Morse chain 
is very exceptional since most of the electronic states ex- 
hibit extended nature despite of the singular continuity 
of the energy spectrum, 5 and it may be interesting to 
investigate the phenomena resulting from these peculiar 
electronic properties. By the way, one of the interesting 
properties of the chains with a singular continuous en- 
ergy spectrum lies in the quantum dynamics of an elec- 
tronic wave packet. The spread of the wave packet in 
the Fibonacci chain is known to be an anomalous power- 
law diffusion; the spread can be diffusive, subdiffusive, or 
superdiffusive, depending on the strength of the system 
parameters. 11 Since the nature of the energy spectrum is 
known to play an important role in quantum dynamics of 
the wave packet, 12 a chain with the singular continuous 
spectrum is expected to exhibit similar dynamic property 
to that of the Fibonacci chain. Then, it may be worth- 
while to test whether the dynamic property of the Thue- 
Morse chain is similar to that of the Fibonacci chain or 
not. We try to answer this question and discuss several 
characteristics distinctive from the case of the Fibonacci 
chain. 

Since Eq. (1) is the Schrodingcr equation for a single 
electron, effects of electron-electron or elcctron-phonon 
interactions, that are inevitable in real systems, arc ne- 
glected. Thus it may be questionable whether the single- 
electron properties obtained from Eq. (1) can be observ- 
able or not when these interactions are present. As a 
step to solve this question, we take into account a kind 
of nonlinear interaction and study the effect of the in- 
teraction on the dynamics of the electronic wave packet. 
We show that degree of dynamic localization is crucially 
dependent on the strength of the hopping energy and the 
interaction parameter. 

This paper is organized as follows: In Sec. II we de- 
duce the trace map of the transfer matrices in the transfer 
model by directly diagonalizing the Hamiltonian matrix 
of the Thue-Morse chain. Then we elucidate the nature 
of the electronic states by calculating the wave functions, 
the Lyapunov exponents, and the Landauer resistivities 
of the states. In Sec. Ill we illustrate the dynamic prop- 
erties of an electronic wave packet by observing the evolu- 
tion of the time-dependent Schrodingcr equation. And, 
in Sec. IV we present the effect of the nonlinear inter- 
action on the dynamics of the wave packet. Finally a 
summary is given in Sec. V. 



II. ELECTRONIC PROPERTIES OF THE 
THUE-MORSE CHAIN 



As a first step to examine the electronic properties in 
the framework of the transfer model, we deduce the trace 
map of the transfer matrices. To do this, we first consider 
the characteristics of the energy eigenvalues obtained by 
directly diagonalizing an N x N Hamiltonian matrix un- 
der a periodic boundary condition. Figure 1 shows a 
plot of E versus I, where / is the order of the chain size 
(N = 2 l ) and T(= t L /t s ) is the strength of the hopping 
energy. From this, we obtain the following characteris- 
tics: (a) The eigenvalues are symmetric with respect to 
E = and there always exists a state with E = 0. (b) All 
the eigenvalues except the two outermost eigenvalues are 
doubly degenerate, and there are (2' _1 + 1) distinctive 
eigenvalues at the Zth-order Thue-Morse chain, (c) For 
I > 2, all the eigenvalues at the rth-order chain remain as 
those at the (I + l)th-order chain, (d) Additional eigen- 
values at the (I + l)th-order chain are located between 
the eigenvalues at the rth-order chain. 

The trace map of the transfer matrices can be deduced 
as follows: First, by noting that property (c) implies the 
fact that (xi-i — 1) is an ingredient of (xi + i — 1), we 
can write xi + i — 1 = f(l)(xi — 1), where xi is the half 
of the trace of M/. Second, by means of properties (c) 
and (d), we draw the eigenvalue tree, as in Fig. 2. Since 
the tree resembles that of the diagonal model, 6 we as- 
sume that f(l) is a power of In order to check this 
assumption, we calculate the energies satisfying 
= to compare them with the energies {Ei +1 } sat- 
isfying xi+\ = 1. Then, we find that becomes a 
half part of {E l+1 }, i.e., {E l+1 } = {E'^Ei}. Thirdly, 
by combining this result with property (b), we find that 
f(l) is a square function of xi-\, i.e, f(l) — 4a; 2 _ 1 . Thus, 
the trace map of the transfer matrices becomes 

x l+1 = 4a; 2 _ 1 (a; ; -l) + l (3) 

with xi = [E 2 - 0| + /2t L t s and x 2 = E 2 [E 2 -2(t 2 L + 
t 2 )]/2t 2 L t 2 +l. 

Using Eq. (3), we can calculate energy eigenval- 
ues for any chain size. For example, let us consider 
the case of T = 0.5. Then, since we can easily ob- 
tain the energies E[ = {±V5}, E' 2 = {±^5 ± \/l7}, 
E 2 = {0,±VT0} and e' 3 = {±[5±[(25±^97)/2] 1 / 2 ] 1 / 2 }, 
we can write the energy eigenvalues up to I — 5 ex- 
actly: E 3 = {0, ± VE, ± VlO}, E 4 = {0, ± 

V5, ± Vl O, ±y / 5±VT7} 1 and E 5 = {0, ± V5, ± 

VW, ± V5 ±yi7, ± [5± [(25 ± V497)/2] 1 / 2 ] 1 / 2 }. 
We have checked that, for higher I, the energy eigenval- 
ues obtained numerically by Eq. (3) are exactly the same 
as those obtained from direct diagonalization method 
within the numerical accuracy. 
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Before going on further, we would like to mention some 
remarks. The first is that Eq. (3) is the same form as that 
of the diagonal model. This is a very surprising result be- 
cause the transfer model has a much more complex pat- 
tern of the transfer matrices compared with those of the 
diagonal model. The second is that we deduced Eq. (3) 
by means of a numerical calculation. Thus, the problem 
of a mathematically rigorous derivation of Eq. (3) still 
remains as an open problem. However, our purpose in 
this paper is not to solve this problem exactly, but to 
elucidate the nature of the electronic energy eigenvalues 
and the eigenstates, let us remain it open. The third 
is that, by using the branching rule, the number of any 
desired subband can be indexed in the same way as in 
Ref. 6. This means that one can also perform a scaling 
analysis on the energy bandwidths as in the case of the 
diagonal model, which is in contrast to the argument of 
Ref. 8. 

From now on, we discuss the nature of the eigenstates 
with energies obtained by Eq. (3). To this end, let us 
first classify the eigenstates into two groups - the degen- 
erate states and the edge states. Since Eq. (3) is exactly 
the same as that of the diagonal model except for the 
initial conditions X\ and x 2 , we affirm that the same be- 
havior as that of the diagonal model also appears in the 
transfer model; degenerate eigenstates at the Ith-order 
chain will exhibit an mth-order-type chain-like extended 
behavior in the wave functions at the (I + m)th-ordcr 
chain. Numerical calculations verify our affirmation. We 
present some examples in Fig. 3. The only exceptions are 
the edge states with energies satisfying xi = 1; the wave 
functions diverge algebraically with increasing the chain 
size, as in the edge states of the periodic chain. 

A little cumbersome calculation enables us to write Mi 

as 



M« = I 



for the degenerate states and 
M t = 



1 + aN bN 
cN 1-aN 



(4) 



(5) 



for the edge states. Here, we have set a — r 2 /4, 
b = T(r 2 /St L )^/2(tl + t 2 s ), c = T(2t L /rt s )b, and r = 
th/ts + ts/tL- Using Eqs. (4) and (5), we can obtain the 
Lyapunov exponent ( defined by 



C= lim Cn = Jim -In||M,| 

AT— >oo N^oo iV 



(6) 



where | | is the modulus of the matrix. Thus, the 
localization length £jv (= Cat 1 ) f° r a chain size N is 
given by £jv = N/\n2 for the degenerate states and 
£jv ~ N/2\nN for the edge states; £/v goes to infinity 
in the infinite limit of the chain size, which indicates the 
absence of exponentially localized states. 

As another test of the nature of the eigenstates, we cal- 
culate the Landauer resistivity. Assuming that a Thuc- 
Morse chain of size N is embedded in the middle of a 



perfect conductor and that an electron with E = 2T cos k 
comes in from the left-hand side of the perfect conductor, 
the dimensionless resistivity can be written as 



p(fc,A0 = ^|t/ 2 i(fc,A0| 2 . 



(7) 



where U2i(k,N) is the (2,1) component of the trans- 
fer matrix relating the coefficients of the incoming and 
reflecting plane waves to that of the transmitted plane 
wave. [See Ref. f3 for the detailed formalism.] Using 
Eqs. (4) and (5), we obtain that \U 2 i(k , N)\ 2 = cot 2 k 
for even I and \U 2 i(k , N)\ 2 = (t s /t L ) 2 cot 2 k for odd I, 
while \U2i(ki, N)\ 2 ~ N 2 . Here, k and k\ are the wave 
numbers that correspond to the degenerate and the edge 
eigenvalues. Thus, we have 



p(k,N) 



1/N 
N 



for 
for 



k = k 
k = k\ 



(8) 



which indicates that the degenerate (edge) states are ex- 
tended (algebraically localized). 



III. DYNAMIC PROPERTIES OF THE 
THUE-MORSE CHAIN 

To study the dynamic properties of the Thue-Morse 
chain, we consider the time-dependent Schrodinger equa- 
tion 



.Olpn 

1 dt 



H tpn 



where 



H ip n = t n+1/n tp n+1 + i„_i,„V> n _i + v n tp n . 



(9) 



(10) 



In integrating Eq. (9), we assume an electronic wave 
packet locates initially at ipn — S n _ no and observe the 
time evolution of the wave packet by employing the 
fourth-order Runge-Kutta algorithm. Besides, we take 
the chain size as N = 16384 (i.e., / = 14) with n = 8192 
and the time step as At = 0.05. Note that, even though 
we use the finite chain size and the fixed boundary condi- 
tion (^o = i>N — 0) , they have no effect on the results of 
the integration since we perform the integration within 
time interval where the wave front does not reach the 
boundary region. We also check the accuracy of the nu- 
merical integration by monitoring the conservation of the 
total probability (i.e., ^ n \ipn\ 2 = 1)- 

In clarifying the dynamic nature of the wave packet, 
the most important quantity is the variance of the wave 
packet 



V(t) = J2(n-nof\Mt)f 



(11) 



which gives a global estimate of the spread of the wave 
packet in space. The asymptotic long-time behavior of 
the variance follows the power-law 
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V(t) ~ t 7 . (12) 

It is well known that 7 = for the localization, < 7 < 1 
for the subdiffusion, 7 = 1 for the ordinary diffusion, 
1 < 7 < 2 for the superdiffusion, and 7 = 2 for the 
ballistic motion, respectively. 

Figures 4 and 5 show the results of V(t) for several val- 
ues of the system parameters T and v (= vl = —vs)- The 
curves in Fig. 4 (5) are obtained within the framework of 
the transfer (diagonal) model. For the curves in Fig. 4, 
we obtain the dynamic index 7 as 7 = 1.942 ± 0.394 
for T = 0.125, 7 = 1.931 ± 0.103 for T = 0.2, 7 = 
1.924 ± 0.039 for T = 0.5, 7 = 1.998 ± 0.001 for T = 1.0, 
7 = 1.853 ± 0.028 for T = 1.25, 7 = 1.924 ± 0.044 for 
T = 2.0, and 7 = 1.840 ± 0.078 for T = 5.0, respec- 
tively. Note that the standard deviation in the calcu- 
lated 7's becomes larger as the value of T deviates from 
T = 1.0, which means that the oscillatory behavior be- 
comes stronger as T — > 00 and T — > 0. However, from 
the trend of the curves, we expect that the oscillatory 
behavior will shrink to zero for sufficiently long times. 
Concerning with the oscillatory behavior, we would like 
to mention a point: Comparing Fig. 4 (5) with Fig. 2 
(6) of Ref. 11, we can see that there is no strong hier- 
archical time evolution of V(t) unlike to the case of the 
Fibonacci chain. This feature may attribute to the weak 
self-similarity in the energy spectrum of the Thue-Morse 
chain. 

Recently, de Brito et al. 14 argued with the diagonal 
model of the Thue-Morse chain that 7 is independent of 
the strength of v and is given by 7 = 1.65. On the con- 
trary, Katsanos et al. 15 argued with the same model that 
7 is strongly dependent on the strength of v and given 
by 7 = 1.5 for large values of v. However, judging from 
our results, both of the arguments are partly correct. As 
can be seen in Figs. 4 and 5, 7 depends on the system 
parameters T and v to a certain extent unlike to the re- 
sult of Ref. 14. However, the dependence seems not to be 
strong unlike to the argument of Ref. 15. Note that this 
is another feature of the Thue-Morse chain distinctive 
from the Fibonacci chain: In the case of the Fibonacci 
chain, 7 decreases continuously from 2 to as the system 
parameters increase, and the diffusion of the wave packet 
is known to be either superdiffusive or subdiffusive, de- 
pending upon whether T _1 < 4 or T^ 1 > 4. 11 However, 
there seems to be no such critical value of T in the Thue- 
Morse chain and the spread of the wave packet seems 
to be always superdiffusive; 1.8 < 7 < 2.0 for T > 1 
and 1.9 < 7 < 2.0 for T < 1 in the transfer model and 
7 ~ 1.50 ~ 1.65 for the diagonal model. The most im- 
portant point to be noted is the fact that the obtained 7's 
are intermediate between the value in the perfect periodic 
chain (7 = 2) and the value in the metallic regime of the 
Harper's model or three-dimensional disordered lattices 
(7 = 1). This feature attributes to the unusual elec- 
tronic properties of the Thue-Morse chain; 7 > 1 reflects 
the extended nature of the electronic states, while 7 < 2 
reflects the singular continuity of the energy spectrum. 



Another important quantity that characterizes the dy- 
namics of the wave packet is the initial-site probability 
defined by R(t) = \i(j no (t)\ 2 , which has been widely used 
in investigating the self-trapping transition of the wave 
packet. 16 However, since the fluctuation in R(t) is too 
large to extract characteristic features of the dynamics, 
we calculate the temporal autocorrelation function 17 

1 /"* 

C(t) = - / R(t')dt, (13) 
* Jo 

which is the time-averaged quantity of the initial-site 
probability. The decay of C(t) for long times is known 
to follow the power-law 

c(t) ~ r s . (14) 

The exponent 8 approaches 1.0 for long times in the pe- 
riodic chain, while it goes to in the disordered chain. 
Figure 6 shows the results of C(t) for several values of 
T in the transfer model, where the power-law decaying 
behavior of C{t) is clearly seen; 5 <~ 0.897 for T = 1.0, 
8 ~ 0.313 for T = 1.25, 8 ~ 0.214 for T = 2.0, and 
8 <~ 0.114 for T = 5.0, respectively. The exponent 8 de- 
creases with increasing T, which implies that the spread 
of wave packet becomes slower as T increases. 
We also calculate the participation number 

^(*)=|eiw*)i 4 } ' ( i5 ) 

which gives a rough estimate of the number of sites where 
the wave packet has a significant amplitude; P(t) = 1 for 
a single-localized state and P(t) — N for a state uni- 
formly extended over N sites. Figure 7 shows the results 
of P{t) for several values of T in the transfer model. P(t) 
increases almost monotonically for T — 1.0, while it in- 
creases with strong oscillation for T > 1.0. The oscil- 
lation in P(t) increases and the values of P(t) rapidly 
decreases with increasing T. The former reflects the sin- 
gular continuity of the energy spectrum and the latter 
implies that it becomes more difficult for a wave packet 
to propagate. 

IV. NONLINEAR INTERACTION AND THE 
DYNAMICS 

To study the effect of the electron-phonon interaction 
on the dynamics of the wave packet we take into account 
a kind of nonlinear interaction 18 and consider the time- 
dependent Schrodinger equation 

*ir = {Ha + Hl)xpn (16) 

where 

#1 = -a|Vn|Vn- (17) 
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Here a gives the strength of the nonlinear interaction. 
Note that the nonlinear interaction we consider describes 
a static short-range electron-phonon interaction, which 
appears in many situations such as, for example, the Hol- 
stein model for polaron theory. 19 

Dynamics of an electronic wave packet is known to be 
sensitively dependent on the nonlinear term in Eq. (16). 
In the case of periodic chains, 20 there occurs self-trapping 
transition when a becomes larger than a certain critical 
value a c . And, in disordered chains, 21 degree of dynamic 
localization weakens with the increase of a. In the mean- 
while, in the case of the Thue-Morse chain, Johansson et 
al. 12 argued with the diagonal model that self-trapping 
occurs for arbitrarily small values of a and that the vari- 
ance of the wave packet increase infinitely with time. 
Thus it may be questionable whether the dynamics of 
the wave packet within the transfer model is different 
from the case of the diagonal model or not. We integrate 
numerically Eq. (16) with setting v = 0. 

Figure 8 shows the variance of the wave packet for 
several values of a. As in the case of the linear case 
(a = 0), the oscillatory behavior shrinks to zero for 
long times. Thus the asymptotic long-time behavior of 
the wave packet is also characterized by a superdiffusivc 
movement in spite of the existence of the nonlinearity, 
which resembles the cases of the linear transfer model and 
the nonlinear diagonal model. 12 Besides, the dynamic in- 
dex 7 is nearly independent of the strength of a, which 
also resembles the case of the nonlinear diagonal model. 
The most interesting point to be noted is that the depen- 
dence of the variance on the strength of a for T > 1 is 
qualitatively different from that for T < 1 : When T > 1, 
the variance decreases with increasing a [Fig. 8(a)], which 
indicates that the nonlinearity resists the propagation of 
the wave packet. In the meanwhile, when T < 1, the 
variance increases with increasing a [Fig. 8(b)], which 
indicates that the nonlinearity assists the propagation of 
the electronic wave packet. Our results show that de- 
gree of dynamic localization of the wave packet in the 
Thue-Morse chain crucially depends on the strength of 
the interaction. 



V. SUMMARY 

In summary, we study the spectral and dynamic prop- 
erties of a one-dimensional Thue-Morse chain mainly 
within the framework of the transfer model. By using 
the method of direct diagonalization of the Hamiltonian 
matrix, we first illustrate the energy spectral properties 
and deduce the trace map of the transfer matrices. The 
trace map is shown to be exactly the same as that in 
the diagonal model. Then we illustrate the nature of the 
electronic states by calculating the wave functions, the 
Lyapunov exponents, and the resistivities of the states. 
All the doubly degenerate eigenstates are shown to be 
extended, while the edge states are algebraically local- 



ized. Next, in the study of the dynamics, we show that 
the asymptotic long-time behavior of an electronic wave 
packet is superdiffusive with the dynamic indices 7 sat- 
isfying 1<7<2;7>1 reflects the extended nature 
of the electronic states, while 7 < 2 reflects the singular 
continuity of the energy spectrum. We also show that the 
dependence of the dynamic index on the system param- 
eters is much weaker than that of the Fibonacci chain 
and that there is no hierarchical time evolution of the 
variance unlike to the Fibonacci chain. Finally, by tak- 
ing into account a kind of nonlinear interaction, we study 
the effect of the electron-phonon interaction on the dy- 
namics of the wave packet and show that degree of the 
localization of the wave packet crucially depends on the 
strength of the interaction. 

In spite of fundamental importance of electron-electron 
interaction, we did not take into account the interac- 
tion in this paper. However, note that, in connection 
with the problem of whether degree of localization of a 
single-electron weakens or strengthens in the presence of 
this interaction, much attention has been devoted to the 
chains with disordered 22 and quasiperiodic 23 potential 
energies. Thus, it is worthwhile to investigate the effect 
of electron-electron interaction on the electronic and dy- 
namic properties of the Thue-Morse chain. This problem 
is on working and will be published elsewhere. 

ACKNOWLEDGMENTS 

This work was financially supported by the Korea Re- 
search Foundation made in the program of 1997-003- 
D00092. 



1 D. Shechtman, I. Blech, D. Gratias, and J. W. Cahn, Phys. 
Rev. Lett. 53, 1951 (1984). 

2 M. Kohmoto, L. P. Kadanoff, and C. Tang, Phys. Rev. 
Lett. 50, 1870 (1983); S. Ostrund, R. Pandit, D. Rand, H. 
J. Schellnhuber, and E. Siggia, ibid. 50, 1873 (1983); Q. 
Niu and F. Nori, ibid. 57, 2057 (1986); J. A. Ashraff and 
R. b. Stinchcombe, Phys. Rev. B 37, 5723 (1988); T. Fuji- 
wara, M. Kohmoto, and T. Tokihiro, ibid. 40, 7413 (1989); 
A. Suto, J. Stat. Phys. 56, 525 (1989). 

3 R. Merlin, K. Bajema, R. Clarke, F. Y. Juang, and P. K. 
Bhattacharya, Phys. Rev. Lett. 55, 1768 (1985). 

4 F. Axel, J. P. Allouche, M. Kleman, M. Mendes-France, 
and J. Jeyriere, J. Physique Coll. 47, C3, 181 (1986); Z. 
Cheng, R. Savit, and R. Merlin, Phys. Rev. B 37, 4375 
(1988). 

5 K. Iguchi, Mod. Phys. Lett. B 1, 29 (1994); A. Bovier and 
J.-M. Chez, J. Phys. A 28, 2313 (1995). 

6 C. S. Ryu, G. Y. Oh, and M. H. Lee, Phys. Rev. B 46, 
5162 (1992); ibid. 48, 132 (1993). 



•5 



7 G. Y. Oh, J. Kor. Phys. Soc. 31, 808 (1997). 

8 A. Chakrabarti, S. N. Karmakar, and R. K. Moitra, Phys. 
Rev. Lett. 74, 1403 (1995). 

9 G. Y. Oh, C. S. Ryu, and M. H. Lee, J. Phys.: Condens. 
Matter 4, 8187 (1992). 

10 J. Bellissard, A .Bovier, and J.-M. Ghez, Commun. Math. 
Phys. 135, 379 (1991). 

11 H. Hiramoto and S. Abe, J. Phys. Soc. Jpn. 57, 230 (1988). 

12 M. Johansson, M. Hornquist, and R. Riklund, Phys. Rev. 
B 52, 231 (1995). 

13 M. Kohmoto, Phys. Rev. B 34, 5043 (1986); Y. Lie and K. 
A. Chao, ibid. 34, 5247 (1986). 

14 P. E. de Brito, C. A. A. da Silva, and H. N. Nazareno, 
Phys. Rev. B 51, 6096 (1995). 

15 D. E. Katsanos, S. N. Evangelou, and S. J. Xiong, Phys. 
Rev. B 51, 895 (1995). 

16 H. Hiramoto and S. Abe, J. Phys. Soc. Jpn. 57, 1365 
(1988). 

17 J. P. Lu, T. Odagaki, and J. L. Birman, Phys. Rev. B 33, 
4809 (1986). 

18 D. Chen, M. I. Molina, and G. D. Tsironis, J. Phys.: Con- 
dens. Matter 5, 8689 (1993). 

19 T. Holstein, Ann. Phys. (NY.) 8, 325 (1959); L. A. Turkc- 
vich and T. D. Holstein, Phys. Rev. B 35, 7474 (1987). 

20 M. I. Molina and G. P. Tsironis, Physica D 65, 267 (1993). 

21 D. L. Shepelyansky, Phys. Rev. Lett. 70, 1787 (1993); M. I. 
Molina and G. P. Tsironis, Phys. Rev. Lett. 73, 464 (1994). 

22 D. L. Shepelyansky, Phys. Rev. Lett. 73, 2607 (1994); S. N. 
Evangelou, S-J. Xiong, and E. N. Economou, Phys. Rev. 
B 54, 8469 (1996); R. A. Romer and M. Schreiber, Phys. 
Rev. Lett. 78, 515 (1997); K. Fraham, A. Miiller-Groeling, 
J-L Pichard, and D. Weinmann, ibid. 78, 4889 (1997); R. 
Romer and M. Schreiber, ibid. 78, 4890 (1997). 

23 H. Hiramoto, J. Phys. Soc. Jpn. 59, 811 (1990); D. L. She- 



pelyansky, Phys. Rev. B 54, 14896 (1996); A. Barelli, J. 
Bellissard, P. Jacquod, and D. L. Shepelyansky, Phys. Rev. 
Lett. 77, 4752 (1996); S. N. Evangelou and D. E. Katsanos, 
Phys. Rev. B 56, 12797 (1997). 



FIGURE CAPTIONS 

Fig. 1. E versus I in the transfer model with T = 0.5. 
Fig. 2. Schematic representation of the eigenvalue tree. 
Fig. 3. \tp n \ versus n for states with (a) E = \/5 

(X! = 0) and (b) E = ^5 + y/Vf (x 2 = 0) with T = 0.5, 
(ipo,ipi) = (0.5, 1), and I = 6. \ip n \ exhibits a third-order- 
type (i.e., LSSLSLLS-type) chain- like behavior for (a) 
and a second-order- type (i.e., L SSL-type) chain- like be- 
havior for (b). 

Fig. 4. log 10 V(t) versus log 10 1 within the transfer 
model: (a) T = 1.0,1.25,2.0,3.0, and 5.0, and (b) T = 
1.0,0.5,0.2, and 0.125 from top to bottom at t = 1000. 
Fig. 5. log 10 V(t) versus log 10 i within the diagonal 
model, where v = 0.1,0.5,1.0,2.0, and 3.0 from top to 
bottom at t = 3000. 

Fig. 6. log 10 C(t) versus log 10 t within the transfer 
model, where T = 1.0,1.25,2.0, and 5.0 from bottom 
to top at t = 1000. 

Fig. 7. P(t) versus t within the transfer model, where 
T = 1.0, 1.25, and 2.0 from top to bottom at t = 200. 
Fig. 8. log 10 V(t) versus log 10 1 within the transfer 
model: (a) T = 5.0 with a = 0, 0.2, and 0.8, and (b) 
T = 0.2 with a = 0.6, 0.4, 0.1, and from top to bottom 
at t = 1000. 



G 




Fig. 1 : Gi-Yeong Oh 



x 2 = 1 



X\ = 


1 








1 T ^ 















i ^2 






1 




I x\ = 


| 














1 rp 2 

, 2 









' x 2 = 


1 


















* rp 2 

1 ^2 













x\ = 






, Xi = 


1 








l . 2 

| tJv 2 



















I x 2 = 


1 








I = 


1 


2 




3 


4 





Fig. 2 : Gi-Yeong Oh 



LSSLSLLS 




Fig. 3(a) : Gi-Yeong Oh 



4 








1 



n 



64 



Fig. 3(b) : Gi-Yeong Oh 




Fig.4 : Gi-Yeong Oh 




Fig.5 : Gi-Yeong Oh 



■10 12 3 4 

logt 

Fig. 6 : Gi-Yeong Oh 




Fig. 7 : Gi-Yeong Oh 



1 r 

(a) 




Fig.8 : Gi-Yeong Oh 



